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Abstract 

We study the correlation between the trigger vr and the associated J/ip on 
near and away sides in Au+Au colhsions at ^/SNN = 200 GeV. In the re- 
gion of trigger momentum pt > A GeV/c, the vr spectrum is composed of 
thermal- shower and shower-shower recombinations in the frame work of the 
recombination model. We consider the azimuthal anisotropy in the quenched 
hard parton distribution and then calculate the elliptic flow parameter V2 of 
charmed mesons {J/ip, and Ds) for different centralities. 

Keywords: Recombination model; di-hadron correlation; elliptic flow 

parameter; charmed meson 
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1. Introduction 

Jet quenching is one of the most informative probes for the formation 
of quark gluon plasma (QGP) in relativistic heavy ion collisions (RHIC). It 
describes that partons produced in the initial collisions with high px lose 
a large fraction of their energies when they propagate through the matter 
before the formation of final state hadrons. The experimental results of the 
suppression of both inclusive hadron [l], Isf and correlated away-side di- 
hadron yields j3] were predicted from this energy loss j5|, loj. Usually the 
energy loss of a jet is caused by two main mechanisms: elastic collisions with 
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deconfined partons and induced gluon radiation. Since the energy loss is 
related to the traverse path length of the jet 0,0], two-particle correlation is 
effective in investigating jet-medium interaction corresponding to the passage 
of a hard or semihard parton |9|]- In Ref . [lol] the correlation of vr-vr in jets 
produced on near and away sides of high pt triggers has been studied in the 
frame work of the recombination model. Now we consider the heavy-flavored 
meson yields and calculate the di-hadron correlation of tt-J /ip in this paper. 

The elliptic flow which is a measure of the anisotropic momentum distri- 
bution of the produced particles is another important observable at RHIC. 
The elliptical flow (^2) results from the initial geometrical anisotropy which 
creates pressure gradients soon after the start of the hydrodynamic evolu- 
tion. Thus, the final V2 is sensitive to the fireball evolution and allows us 
to probe the dynamical properties of the dense matter produced in the col- 
lisions. Recently PHENIX Collaboration has measured the J /if) elliptic flow 
parameter V2 in Au-|-Au collisions at ^Jsnn = 200 GeV In this paper, 
we calculate f 2 in the recombination model and compare our results with the 
experimental data. 

This paper is organized as follows. We calculate the di-hadron correlation 
of TT-J/ip and V2 of charmed mesons in Section 2 and 3, respectively. And a 
brief summary is given in the final section. 



2. Di-hadron correlation of tt-J/i/) 

Based on the recombination model we have calculated the transverse 



momentum spectrum of J/ip which fits the experimental data well [12|, and 
also have predicted pt spectra of other heavy flavored mesons for different 
centralities [HI]. In this section we apply the shower parton distributions 
(SPD) to calculate the di-hadron correlation of n-J/tp on near and away 
sides. We use p as a generic symbol for the transverse momentum of any 
hadron, pt for the trigger particle (vr), pa for the associated particle {J/ip) on 
the near side and pb for the associated particle (J/ip) on the away side. 

In the recombination model the meson production is expressed as the sum 
of three terms, TT (pure thermal), TS (thermal-shower) and SS (shower- 
shower), since there are two components of parton sources: thermal (T) 
partons and shower (S) partons which originate from hard partons. Thus we 
have 
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The shower-shower recombination term {SS contribution) is equivalent 
to the fragmentation function (FF) [lo| 
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is the distribution of hard parton i in Au+Au collisions at ^JsNN = 200 
GeV after traversing the medium of length L. The distribution fi{k) = 
dN^'^'^'^ / (Pkdy of hard parton i just after hard scattering at mid-rapidity can 



be found in Refs. |14|| and [15|. /3L is the explicit dynamical medium factor 
to describe the energy loss effect and its value for different centralities is 
determined by fitting the single-pion inclusive distribution 10|. Df^ is the 



FF of parton i splitting into meson M and the results of FFs for vr and J /il) 



are obtained in Ref . [16j] and Refs. [l7|, ll8| , respectively. 



The thermal parton transverse momentum distribution density is given 
in Ref.[l9j 
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It's assumed that hadronization occurs at r = 5 fm with temperature T = 
175 MeV in the parton phase 19|], which is consistent with predictions of the 
phase transition temperature at vanishing baryon chemical potential from 
lattice QCD 20|. At = PnT^ is the transverse area of the parton system with 
the radius po = 9 fm [19|. Using the parameterized SPD Sj{z) [l^, |2l|, we 
can determine the distribution of shower parton j with transverse momentum 
qi in central Au-|-Au collisions as [iq| 

dq 

q 
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where 



1 rqe'^^ ^h. 



(6) 



is a little different from Eq.(l3]) [12] 
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Then the TS term can be calculated as 
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where 



with the momenta of the two constituent quarks qi = xp and g2 = (1 — x)p. 
4>m{x) in Eq.(IH]) is the wave function of the meson in the momentum space. 
The wave functions for tt and J/^ are determined in Refs. 2^ and (l2| . 
respectively. 7a and stand for the fugacities of the constituent quarks a 
and b. The fugacities of hght quarks are 7« = 7d = 1, 7s = 0.8 [l9| and 
for charm quark 7c = 0.26 which is obtained by fitting the J/ip transverse 
momentum spectrum 12j. Since every quark has 3 color and 2 spin degrees 
of freedom, we use the meson degeneracy factor Cm = (3 x 2)^. 

We consider the di-hadron correlation at intermediate trigger momentum 
{pt > 4 GeV), thus approximately the TT contribution to the trigger and jet 
can be neglected. Then the correlation of the particle J/ ip with pa associated 
with the trigger particle vr with pt on the near side is 



PtdptPadpa 



Y: [ ^F,m[fs\q,p,) + ^Dn^-^)]Ts'^\q-p„p^) 
i J (1 PH q 

+fs\q-p^,p,)^Dil^{^) + ^)}, (9) 
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where the di-hadron FF (-21,-^2) is assumed as [10 

A(^i,^2) = l[D-{z,)D'/^{-^) + D^{-^)D'/^{z,)]. (10) 

Z i — Zi i — Z2 

The di-hadron distribution of the tri gge r vr and the associated J/ip with pb 
on the away side is calculated in Ref . [10| as 

x|f§"*(,',ft) + 4,B;"*(^)l (11) 

Pbq q 

with qQ = Max(ge~^^,pb). Here, q' is the momentum of the recoil parton in 
the trigger jet. 
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Then the near-side (away-side) yield per trigger for the trigger momentum 
in a narrow range Apf (Ap^ — > 0) around pt is 



near(away)^/ dp.^^f-^ / dp,^ (12) 

'^^ iApt Pa{Pb)dptdpa{pb)' J Apt dpt 

where dN^^/dpt is the trigger pion distribution that excludes the TT compo- 
nent of the inclusive distribution. 

We use c to denote the centrality, for example, c = 0.05 stands for 0-10% 
centrality. The energy loss factor /3L is dependent on c and its values can be 



found in Ref.[10[- With the known parameters of /3L we can calculate the 
yields for different centralities. The results of the yield of J/ ip on near and 
away sides per trigger with three values of trigger vr transverse momentum 
Pt in Au+Au collisions for 0-10% and 40-50% centralities at ^/sNN = 200 
GeV/c are shown in FigJTJand Figj2l respectively. For the central collisions 
(c = 0.05), the near-side yield increases with pt for pa > 3.8 GeV/c, while 
this trend is reversed in the region of low Pa- The away-side yield becomes 
higher as pt increases. We find that the shape of the Pa distribution does not 
change with the centrality, while the trend of the away-side yield becomes 
smoother and similar to that of the near-side in more peripheral collisions. 

In order to clarify the dependence of the yield on the centrality, the near- 
side yield versus c is shown in Fig JS] for pt = 4 and 6 GeV/c and pa = 2 
and 6 GeV/c. The result is nearly a constant in c, which has been testified 
in Ref.[loj. The calculation of the average distance (t) the parton travels to 
reach the near-side surface [lo| suggests that the hard scattering point is in a 
layer roughly 13% of the medium size L inside the surface and is insensitive 
to L. On average 15% of the parton energy is lost to the medium, which is 
also independent on L. Thus the centrality dependence of the near-side yield 
is negligible. The away-side yield increases when c is raised from 0.05 to 0.75 



for pb > 3 GeV/c as shown in FigJH It has been explained in Ref. [10| when a 
particle on the away side is required, the scattering point cannot to be too far 
from the surface of the away side and the point is pulled closer to the surface 
with higher c. So when the nuclear overlap is smaller, it is easier for the 
recoil jet to reach the away side to produce a particle. The calculation of the 



ratio of recoil parton average momentum {q') / {k') in Ref. [10| also suggests 
that the energy loss of the recoil parton becomes smaller with increasing c. 
That means in higher c collisions the recoil parton emerges at the away-side 
surface with a larger momentum q' at given initial momentum k'. Then, at 



5 



low pb this leads to a visible smaller momentum fraction x = pb/q' in Eq.f llip . 
The FFs shown in Ref.|2l| reflect that the probability of J/ip produced at 
low X is very small, which results in the lower yield accordingly. Thus the 
behavior of the dependence on the centrality is quite different in the region 
of low Pb- 



The difference between tt-J /ip and tt-tt correlation calculated in Ref.(lO 



caused by the significantly different FFs for J /if) and tt. One can take a clear 
comparison of the FFs in Refs.|12j] and [2l[. It's necessary to note that in 
Ref. [lO| the shower partons can be initiated by hard partons i = u,d, s, u, d, s 



and gluon. Thus in vr-vr calculation the FFs for pion come from three terms 
D^, D1 and D'^ where v denote the valence quark, s the sea quark and g 
the gluon. In tt-J /ip correlation most contributions are from terms in which 
i can be c and g since the FFs of light quarks splitting into J / are assumed 
to be zero in Ref. 17|. The FF D^^^ increases up to the momentum fraction 
X ~ 0.75 [x ~ 0.25 for Dg^^) and then decreases, while FFs for tt decrease 
monotonically with x. It means that the probability of a parton splitting 
into TT is much higher at low x but for J/ ip the larger probability takes place 
at intermediate x. This contrary trend of FFs at low x results in the different 
behavior of the near-side yield between tt-tt and n-J/ip for pa < 3.8 GeV/c. 
This difference also affects the away-side yield at low pb when c increases 
with smaller PL, which is shown in the lower panel of Figj2l The yield for 
TT-J/ip is 4 ~ 5 orders of magnitude lower than that of tt-tt since the FFs for 
J/ip are much smaller than those for tt and the hard parton distribution of c 
is about 3 orders lower than those of light quarks. 

What we have calculated is the transverse momentum spectra dN/pdp 
of TT and J/ip averaged over all azimuthal angle 0. Because of the average 
over (f) in the spectra, only the effective traversing length /3L appears in 
the calculations for the medium effect. In the next section, we consider 
the azimuthal anisotropy in the spectra of charmed mesons and discuss the 
elliptic flow parameter. 



3. V2 of charmed mesons 

We discuss the azimuthal anisotropy of charmed meson momentum distri- 
bution by considering the elliptic flow V2- We assume that the distribution of 
thermal parton is isotropic so that the recombination of TT is independent 
of azimuthal angle (0). Then the dependence comes from the degraded 
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hard parton distribution, which is caused by the dependence of the travers- 
ing length in the medium. Both TS and SS components are proportional 
to the distribution of the hard partons. Thus we can rewrite the meson 
production (Eq.([T])) as 

c 



pdpdcj) 



J2 [ ^F,i<l,M[fS{g,p) + ^DriP)] + ^, (13) 
J q p^q q pdp 



where Fi{q,(f),c) is the probability of a hard parton i with momentum q at 
azimuthal angle in Au+Au collisions with centrality c at ^/snn = 200 
GeV. 

For a hard parton created with the distribution fi{k) at the creation 
point, the initial momentum k changes into q after traversing an absorptive 
distance ^. The corresponding distribution is expressed by the momentum 
degradation factor G{k, q, ^) 

F^{q,0 = j dkkf,{k)G{k,q,0. (14) 

In terms of ^, G(k, q,^) can be written as a simple exponential form G{k, q,^) = 
qS{q — ke~^) [22[. If P(^, 0, c) is the probability of having a dynamical path 
length ^ for a parton directed at 0, the hard parton distribution is obtained 
after carrying out the integration over ^ 

F,(g, 0, c) = f d^P{^, 0, c)F,(g, 0- (15) 



In Ref. [22|, a scaling behavior of P(^, 0, c) is found for the dependencies on 
and c for pion production. P(^, 0, c) can be written as a universal function 
in terms of a scaling variable z = ^/^(0, c) 

p(e,0,c) = ^(z)/e(0,c), (16) 

where ^ is the mean dynamical length. The scaling function is parameterized 
by 

iIj{z) = C (1 - Cr/B{ar + 1, + 1), C = z/2A (17) 

with ai = 0.37 and 02 = 0.81, where B{a,P) is the Beta function. We can 
get the results of the average dynamical path length ^ as a function of and 



c in Ref. 22 . 
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The elliptic flow is quantified by the second Fourier coefficient V2 of the 
meson azimuthal distribution 

■ = A{pT)[l + 2v2ipT)cosi2(l))]. (18) 



PTdpTdcp 

Then we can get the elliptic flow parameter V2 of charmed mesons {J/ip, -D°, 
Ds) dependent on transverse momentum for different centralities. 

The results of J ftp V2 for different centralities are exhibited in Figj5]with 



the experimental data from PHENIX Collaboration The calculated 

results are in agreement with the data within errors except the negative 
point at pt = 1-5 GeV/c. V2 increases monotonically up to — 4 GeV/c 
and then starts to saturate to constants 0.04 — 0.14 from c = 0.05 — 0.45, 
depending on the centrality. And the values of f 2 get larger from 0.05 — 0.45 
centralities at given px. The trend of pt and centrality dependencies of w 2 is 



similar to that of other charge hadrons shown by PHENIX in Ref . |23 

We also predict V2 of and Dg and the results for are shown in 
FigE] and the comparison of V2 for the three different charmed mesons at 
c = 0.05 is exhibited in the lower panel. Obviously, V2 of or Ds is much 
larger than that of J/ip. This feature results from the different constituent 
quark mass in the three mesons. The dependence of momentum distribution 
on (f) in Eq. flT^ comes from two terms TS and SS with the hard parton 
distribution Fi{q,(l),c). The thermal parton distribution of light quarks {u, 
d and s) is much higher than that of heavy quark [12[. Thus, for D*^ the 
main contribution to TS comes from TuSc- And its transverse momentum 
is mainly from that of the charm quark. Roughly one can take p^o ~ p^.. 
For J /ip^ since its two constituent quarks have the same mass, the transverse 
momentum for both c and c in J are the same as half of that of J /i). Thus 
one would expect that V2^{pt/2) ~ And V2^^{pt/2) is shown in 

the lower panel of Figl6] by the dashed line with symbols, f 2 for and Dg 
is almost the same and larger than that of J/ip at the same pt- 

The number of quark (NQ) scaling of V2, V2/n versus Pr/n, has been 
observed at RHIC 24, 25||. The scaling lends strong support to the finding 



that the collectivity develops in the partonic stage at RHIC. Quark coales- 
cence model has been used to explain the fiavor dependence of V2 [ioj. In this 
model hadrons are combined by the coalescing quarks and then f 2 of hadrons 
can be expressed in terms of f 2 of quarks, where there are no difference in 
collectivity among u, d and s quarks. Such a scaling implies that the mo- 
mentum of the hadron is simply the sum of the momentum of the coalescing 
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quarks. In our model the momentum of or Dg comes mainly form charm 
quark due to the unequal constituent quark masses ^ muinis). Thus 
we get different V2 for J/ip and or Ds- Our prediction can be tested by 
further experiments. If the experimental data get the larger value of v 2 for 
as we have expected, it offers a good proof for the success of the recom- 
bination model. So the violation of the NQ scaling also cannot prove the 
non-formation of QGP in the collisions. 

4. Conclusions 

The correlation of -K-J/ip and the elliptic flow parameter V2 of charmed 
mesons are studied in the frame work of the recombination model. In the 
region of trigger particle momentum > 4 GeV/c, the contribution to the 
meson production comes from TS and TS recombination (TT contribution 
can be neglected). The away-side yield per trigger increases with the trigger 
momentum pt. And the difference between tt-tt and t^-J /ip correlation takes 
place in the near-side yield in the region of Pa < 3.8 GeV/c and the away- 
side yield in collisions with higher centrality at low Ph, which is caused by 
the different behavior of the fragmentation functions for tt and J/ip. 

The elliptic flow parameter V2 of charmed mesons for different centrahties 
is obtained by considering the azimuthal anisotropy in the hard parton distri- 
bution Fi{q, 0, c). V2 increases up to pt ~ 3.5 — 4 GeV/c and then saturates. 
V2 of and Dg is almost the same in the whole pt region and larger than 
that of J/t/j. 
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Figure 1: Distribution of associated J/%1> (pa) on the near side of a jet triggered by a pion 
with three momenta (pt) for c = 0.05 and c = 0.45. 
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Figure 2: Distribution of associated J /ip (pf,) on the away side of a jet triggered by a pion 
with three momenta {pt) for c = 0.05 and c = 0.45. 
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Figure 3: Yield per trigger in the near-side jet as functions of centrality c with Pa and pt 
in GeV/c. 
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Figure 4: Distribution of associated J/ip {pb) on the away side of a jet triggered by a pion 
for three centraUties at = 6 GeV/c. 
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Figure 5: V2 for J/tp as a function of px in Au+Au collisions at ^/sN7f = 200 GeV for 
different centralities. The data shown for three pr^s are from 11 [. 
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Figure 6: V2 for as a function of pT in Au+Au collisions at ^smn = 200 GcV for 
different centralities and the comparison of the elliptic flow parameter V2 for three charmed 
mesons (c=0.05). 
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